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Abstract 


Human coronaviruses (HCoV) are RNA viruses that cause respiratory tract infections with viral 
replication of limited duration. The host and viral population heterogeneity could influence clinical 
phenotypes. Employing long RT-PCR with IIlumina sequencing, we quantified the gene mutation load 
at 0.5% mutation frequency for the 4,529 bp-domain spanning the Spike gene (4,086 bp) of HCoV- 
OC43 in four upper respiratory clinical specimens obtained during acute illness. There were a total of 
121 mutations for all four HCoV samples with the average number of mutations at 30.3 + 10.2, which 
is significantly higher than that expected from the Illumina sequencing error rate. There were two 
mutation peaks, one at the 5’ end and the other near position 1550 in the $1 subunit. Two coronavirus 
samples were genotype B and two were genotype D, clustering with HCoV-OC43 strain AY391777 in 
neighbor — joining tree phylogenetic analysis. Nonsynonymous mutations were 76.1 + 14% of mutation 
load. Although lower than other RNA viruses such as hepatitis C virus, HCoV-OC43 did exhibit quasi- 
species. The rate of nonsynonymous mutations was higher in the HCoV-OC43 isolates than in hepatitis 
C virus genotype la isolates analyzed for comparison in this study. These characteristics of HCoV- 
OC43 may affect viral replication dynamics, receptor binding, antigenicity, evolution, transmission, 


and clinical illness. This article is protected by copyright. All rights reserved 
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1. Introduction 


Human coronavirus-OC43 (HCoV-OC43) is within the Betacoronavirus genus (family Coronaviridae) 
and is an enveloped, positive sense, single-stranded RNA virus [Belouzard et al., 2012; Bosch et al., 
2003]. There are five known genotypes and other betacoronaviruses include severe acute respiratory 
syndrome (SARS) CoV and Middle East respiratory syndrome (MERS) CoV [Ren et al., 2015]. HCoV- 
OC43 is prevalent among humans and genotype D has been prominent in recent years [Lau et al., 2011; 
Ren et al., 2015; Zhang et al., 2015]. Little is known about how HCoV-OC43 genotypes persist in 
human populations, but continuous adaptation by viral antigenic genes in the Spike protein through 
genetic drift may be necessary. The Spike protein is the major antigenic protein and is under selection 
pressure by the host immune response; it is important for host range and tissue tropism. It is cleaved 
into S1 and S2 subunits for receptor binding and membrane fusion. The N-terminal domain of the S1 
subunit is responsible for sugar receptor binding and the $2 subunit is responsible for fusion of viral 
and host membranes [Li, 2016]. The S1 subunit is more divergent in sequence and the S2 subunit is 


more conserved [Belouzard et al., 2012; Masters and Perlman 2013; Ren, 2015]. 


Human coronaviruses cause the common cold and influenza-like illnesses, but can be associated with 
more severe illnesses such as pneumonia, exacerbations of asthma and chronic obstructive pulmonary 
disease, croup and bronchiolitis. In patients with chronic obstructive pulmonary disease studied during 
the 1998 to 1999 influenza season, 13.5% of illnesses were associated with HCoV-229E and HCoV- 
OC43 infection, while in another study between 2009 and 2013, 19% of acute respiratory illnesses in 
patients with cardiopulmonary diseases and 21.5% in healthy young adults were associated with HCoV 
[Gorse et al., 2003, 2006, 2009, 2015]. Coronavirus-associated illness was less severe than influenza 
but was associated with multiple respiratory and systemic symptoms, and hospitalization [Gorse et al., 


2009]. HCoV-229E and HCoV-OC43 infection rates of 2.8% to 26% in healthy young and elderly 
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adults, high-risk adults, and hospitalized patients were reported during the winters of 1999 to 2003 and 


they contributed to medical disease burden [Walsh et al., 2013]. 


Little is known about the degree of heterogeneity of HCoV-OC43 viral quasi-species present in upper 
respiratory secretions. If present, this may help explain persistent incidence of HCoV-OC43 infections 
in human populations, if the mutational changes result in antigenic drift. This might allow escape from 


host immunity and contribute to virus infectivity and pathogenicity. 


In the current study, we combined RT-PCR and Illumina sequencing to measure the diversity of 
HCoV-OC43 Spike gene quasi-species through direct count of the Spike gene mutations, determination 
of percent nonsynonymous mutation rates and comparison of these rates to hepatitis C virus (HCV), 
which is in the genus Hepacivirus, family Flaviviridae. HCV is an RNA virus with a heterogeneous 


population of quasi-species in chronically infected patients [Fan et al., 2009, 2010; Wang et al., 2014]. 


2. Material and Methods 


2.1 Patient Samples 


We studied nasal and oropharyngeal swab specimens that were obtained from each of four patients 
early during symptomatic acute respiratory illness and positive for HCoV-OC43 nucleic acids by 
multiplex RT-PCR [Gorse et al., 2015]. Serum and nasal wash specimens were obtained at the time of 
acute illness and 3 to 4 weeks after illness onset. They were assayed by enzyme-linked immunosorbent 
assay for serum IgG and nasal wash IgA antibodies to tissue culture-adapted HCoV-OC43 (American 
Type Culture Collection #VR-1558, GenBank: NC_005147.1) that was inactivated by psoralen 
compound and long-wavelength ultraviolet light, as described [Gorse, 2009, 2010, 2015]. Severity of 
acute respiratory illness was measured by two scores: a self-reported visual analogue scale of overall 
illness severity, ranging from 1 (mildest) to 10 (most severe), and a severity of influenza-like 


symptoms and signs score that was the sum of 16 symptoms and signs that were graded on a scale of 0 
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(absent) to 15 (most severe) with a maximum score of 240, as described [Arden et al., 1988; Gorse et 
al., 2009, 2015]. Respiratory and systemic symptoms of the acute illness were recorded. The patients 
gave written informed consent and the study was approved by the Institutional Review Boards at the 
VA St. Louis Health Care System and Saint Louis University. Two recombinant clones from a previous 
study, #1701 and #1709, each containing a 9,022 bp HCV insert, were used to estimate potential errors 
associated with Illumina sequencing [Fan and DiBisceglie, 2010; Wang et al., 2014]. Also, 19 HCV 
genotype la samples from an earlier report [Ren et al., 2015] were available for re-analysis and 


comparison in the current study. 


2.2 RNA Extraction, RT-PCR and Illumina Sequencing 


Total RNA from each nasal and oropharyngeal swab specimen sample was purified using the QlAamp 
Ultrasens Virus Kit (Qiagen, Valencia, CA) according to the manufacturer’s procedures. RT-PCR was 
then applied to amplify a 4,529 bp amplicon spanning the full-length spike gene (4,086 bp). In brief, 
10.6 uL of extracted RNA was mixed with 9.4 uL RT matrix consisting of 1x SuperScript III buffer, 10 
mM DTT, | uM OC43RI1 (reverse primer, 5'-TGC CCC ACA TAC CAC ACA G-3’, position 28164- 
28182, numbering is according to HCoV-OC43 strain, GenBank accession number: AY391777), 2mM 
dNTPs, 20 U of RNase OUT recombinant Ribonuclease Inhibitor, and 200 U of SuperScript II 
Reverse transcriptase (Life Technologies). After 75-min. incubation at 50°C and subsequent 
inactivation, an aliquot of 5uL of RT reaction was applied for the first round of PCR that contained 1x 
GC enhancer (New England Biolabs), 1x Q5 buffer (New England Biolabs), 1.6 mM dNTPs, 0.4 uM 
OC43F1 (forward primer, 5'-GTA CAG GTT GTT GAT TCG CG-3’, position 23210-23229), 0.4 uM 
OC43RI1 and 1.6 U Q5 High Fidelity DNA Polymerase (New England Biolabs). After initial heating at 
94°C for 1 min., cycle parameters were programmed as the first 10 cycles of 94°C for 30 sec., 65°C for 
30 sec. and 68°C for 5 min. followed by 20 cycles of 94°C for 30 sec., 60°C for 30 sec. and 68°C for 5 


min. with a 2 sec. autoextension at each cycle. Two uL of the first round of PCR product was used for 
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the second round amplification with primers OC43F2 (forward primer, 5’-TCT GGC CTC TCT ACC 
CCT ATG GC-3’, position 23439-23461) and OC43R2 (reverse primer, 5'-CTT GAT TAC GGC ACC 
AAG CAT GAC-3’, position 27944-27967), under the same cycle parameters as the first round of PCR. 
Product at expected size was gel-purified using QIAquick PCR purification Kit (Qiagen) and 
quantitated. About 4 to 5 ug of purified DNA product was subjected to library construction. The 
fragment library was constructed using Illumina Nextera XT DNA library preparation kit, and followed 


by Illumina sequencing on NextSeq 500 machine with 1 x 250 bp read output. 


2.3 Sequence Data Analysis 


We first estimated the error rate associated with Illumina sequencing using two recombinant HCV 
clones. In doing so, raw sequence reads in fastq format were first filtered in PRINSEQ (v0.19.5) for 
quality control, including read length >70 bp, mean read quality score >25, low complexity with DUST 
score <7, ambiguous bases <1% and all duplicates [Schmieder and Edwards, 2011]. Filtered reads were 
mapped onto HCV genotype la prototype strain H77 (GenBank accession number AY009606) using a 
gapped aligner Bowtie 2 [Langmead and Salzberg, 2012]. Mapped files were then converted into 
binary format (BAM), sorted and indexed in SAMtools [Li et al., 2009] followed by local realignment 
and base quality recalibration in Genome Analysis Toolkit (GATK) [DePristo et al., 2011]. Next, by 
converting post-alignment BAM files into mpileup format in SAMtools, the consensus sequence for 
each clone was called in VarScan (v2.2.3) with the settings of >1,000x coverage, =25 base quality at a 
position to count a read and >50% mutation frequency [Koboldt et al., 2012; Quinlan and Hall, 2010]. 
The entire pipeline was repeated using individual consensus sequences. Mutations were called at each 
position in VarScan under the setting of 0.5% frequency and base quality from 15 to 40, followed by 


manual check in the Integrative Genomics Viewer [Koboldt, 2012]. 


Using the value of base quality to define a mutation from above analysis, similar procedures were 


applied to four patient samples. The HCoV-OC43 strain (GenBank AY391777) was used as the 
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reference at initial mapping. Over the entire coronavirus Spike gene, the mutation load, the total 
number of mutations at a given site, was counted through sliding windows, size = 300 bp, overlap = 
100 bp. Finally, under the frame of full-length HCoV Spike gene (4,086 bp), the nature of each 
mutation, either synonymous or nonsynonymous, was determined using a custom script [Van 


Belleghem et al., 2012]. 


2.4 Phylogenetic Analysis 


The consensus full-length HCoV spike sequences from four patients and reference sequences retrieved 
from GenBank were used for phylogenetic analysis. The tree was constructed using neighbor-joining 
approach under nucleotide substitution model of maximum composite likelihood in MEGA program 


(version 5.2) [Tamura et al., 2011]. 


2.5 Statistical Analysis 


Statistical analyses were done with either two-tailed, unpaired Students test or Chi-square. When 
applicable, data were expressed as mean value and standard deviation. P<0.05 was considered 


statistically significant. 


2.6 Data Availability 


Raw sequence data in fastq format from all four patient samples were archived in NCBI Sequence Read 


Archive (SRA) under SRA accession number SRP071020. 


3. Results 


3.1 Clinical Characteristics of HCoV Infections and Antibody Responses. 


Samples 3 and 4, both genotype D, were collected within a month of each other in December 2010 and 


January 2011 from two older patients with significant acute respiratory and systemic symptoms (Table 
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1). The patients had underlying chronic cardiopulmonary diseases and diabetes mellitus. The two 
illnesses were associated with greater than a four-fold increase in nasal wash IgA antibody titers but 
only one with at least a four-fold increase in serum IgG antibody titer to HCoV-OC43, comparing acute 


illness to convalescent specimens collected 3 to 4 weeks after illness onset (Table 1). 


Samples from subjects 2 and 6, both genotype B, were collected about two years apart in January 2010 
and March 2012 from a younger patient without underlying chronic illnesses and an older patient with 
cardiac disease and diabetes mellitus. Both had acute respiratory and systemic symptoms that may have 
been less severe than those reported by the two patients with genotype D isolates (Table 1). One of the 
two illnesses with genotype B viruses was associated with a greater than four-fold increase in nasal 
wash IgA antibody titer to HCoV-OC43, but neither had a four-fold rise in serum IgG antibody titer to 
HCoV-OC43, comparing acute illness to convalescent specimens collected 3-4 weeks after illness 


onset (Table 1). 


3.2 Quantitation of HCoV-OC43 Mutation Load 


The raw data output indicated 70.1% of bases read had a quality score greater than 30. Interpretation of 
the distribution statistics of base quality scores over read length resulted in trimming the read length at 
the 3’ end by 6 to 10%. The final results of the quality control are shown in Supplemental Table. The 
large output gave a very deep base coverage for each HCoV sample, the average was 94,899 + 21,405 


(Supplemental Fig. 1). 


To estimate the error rate associated with library construction and Illumina sequencing, mutations were 
called from two recombinant HCV clones, #1701 and #1709, under a range of base quality settings. 
Even if a mutation was counted under the base quality as low as 15, there were no differences in the 
consensus sequences derived either from Illumina or from gene-walking Sanger sequencing (data not 


shown). However, the number of individual mutations had a sharp drop from the base quality 25 to 30 
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(Supplemental Fig. 2). Under the conditions of 0.5% mutation frequency and base quality score of =30, 
there were a total of 31 mutations in the two HCV clone samples, suggesting an error rate of about 1.76 


mutations per kb. 


Applying the same criteria for the coronavirus samples, a total of 121 mutations for all four samples 
were identified with an average number of mutations of 30.3 + 10.2 (range: 20 to 40 mutations per 
sample), which is significantly higher than that expected from the Illumina sequencing error rate (121 
vs. 28.76 mutations, P = 4.2 x 10°'*). Nonsynonymous mutations accounted for between 61% and 90% 
of the total mutations (Table 2). No deletions or recombinations were detected. Of the 121 viral 
mutations, those with frequencies greater than 2% occurred at six positions: Spike gene position 


numbers 79, 81, 1229, 1859, 2244 and 2858 (Table 3). 


Using a sliding window analysis with window size = 300 bp and overlaps = 100 bp, there were two 
mutation peaks, one at the 5’ end and the other at about Spike gene position number 1550 (Figure 1). 
Using consensus sequences, a Neighbor-joining tree showed phylogenetically that two of the subjects’ 
samples were genotype D and two were genotype B strains of HCoV-OC43 (Figure 2 A and Figure 2 
B). The two genotype B samples both had high mutation frequencies at spike gene nucleotide positions 
79 (36.92% and 45.12%) and 81 (7.96% and 7.07%), whereas the two genotype D samples had 


mutation frequencies of 0.92% and 0.82% at position 81, and less than 0.5% for position 79. 
3.3 Comparison of HCoV and HCV. 


A previously reported study found that HCV genotype 1a patients experiencing relapse after antiviral 
treatment (n=19) had a higher average total mutation load measured through 454 sequencing [Ren et 
al., 2016]. These samples were re-analyzed for the current study using a base quality score > 30 rather 
than 25 in the earlier report [Ren et al., 2016]. The average mutation load in HCV patient isolates was 


significantly higher than in HCoV-OC43 patient isolates (296.2 + 102.2 vs. 30.3 + 10.2, P=7.7 x 10°). 
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However, nonsynonymous mutations as a percentage of the total mutations were higher among the 
HCoV-OC43 isolates than among the 19 HCV genotype 1a patient isolates (76.7 + 14% vs. 26 + 8%, P 


=35% 107), 
3.4 Consensus Amino Acid Sequences for HCoV-OC43. 


The alignment of Spike genes of consensus amino acid sequences show differences between the four 
clinical strains and the prototype AY391777 strain particularly at the N and C terminal ends of the S1 
subunit, and at the S1 and S2 subunit cleavage site, (a.a. 762-766) with a smaller number of amino acid 
differences in the S2 subunit (Figure 3). The high genetic mutation rates at sites 79 and 81 correspond 
to several amino acid changes compared to the prototype strain between amino acid positions 20 and 


31. 
4. Discussion 


Through Illumina sequencing, we have introduced high-resolution HCoV-OC43 Spike gene mutational 
load to quantify viral quasi-species population diversity. The experimental method allows the role of 
HCoV-OC43 Spike gene heterogeneity to be investigated in a manner not previously reported. The 
HCoV-OC43 Spike gene does have quasi-species, but the magnitude is almost 10 times lower than the 
quasi-species found in HCV, while the nonsynonymous mutations take a higher percentage of the total 
mutation load for HCoV-OC43 than HCV genotype la. It should be noted that the Spike gene was 
amplified with the gene-specific primers. Although these primers are located in the conserved HCoV 
HE and NS2 domains, the missing of potentially heterogeneous viral variants during the amplification 
cannot be excluded. As a consequence, quasi-species diversity of HCoV may be underestimated in the 
current study in comparison to the use of degenerate primers or primer-independent approaches. Viral 
replicative dynamics, population size, and host immune responses may contribute to this observation, 


and the finding has implications in terms of HCoV evolution and treatment. High mutation load 
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facilitates the evolution of viral populations in response to external pressure. Given their slightly 
deleterious nature [Eyre-Walker et al., 2002]; however, excessive nonsynonymous mutations can be 
detrimental to such adaptive evolution [Peck, 1994]. Consequently, accumulation of deleterious 
mutation load during viral infection may contribute to self-limited active infection. In the case of 
HCoV-OC43, this could contribute to short, self-limited respiratory tract infections compared to HCV, 
which is a chronic, systemic viral infection with ongoing viral replication. However, genetic diversity 
manifested by a cloud of quasi-species in some viral systems can be linked to pathogenicity and shown 
to allow adaptation to new environments such as infection of new hosts [Vignuzzi et al., 2006]. If 
clonal virus populations are present in human infections with Middle East Respiratory Syndrome 
(MERS)-CoV [Briese et al., 2014; Cotten et. al., 2013], this could increase the chances of epidemic 
spread of a particularly virulent strain [Gardner and McIntyre, 2014]. Infections with a cloud of viral 
strains may enhance spread, for instance, by increasing the chance of a virulent strain being 
transmitted, or reduce it due to defective and less virulent strains being present. MERS-CoV genome 
sequences from dromedary camels indicate presence of quasi-species in single samples [Briese et al., 
2014]. Recent deep sequencing analyses reported intra-patient viral heterogeneity during a 2015 
outbreak of MERS-CoV and the possibility that host immune response provided selection pressure to 
favor genetic heterogeneity [Park et al., 2016; Kim et al., 2016]. Intra-patient viral heterogeneity may 


have contributed to transmissibility [Park et al., 2016]. 


Our four clinical isolates were of the genotypes B and D which have been described as circulating in 
recent years rather than genotype A [Lau et al., 2011]. Alignment of spike gene consensus amino acid 
sequences compared to the AY391777 prototype strain sequence (genotype A) indicates the majority of 
changes were in the S1 subunit. The proteolytic cleavage site of S1 and S2 subunits for the four clinical 


strains in our study had the RRSRR motif rather than the RRSRG motif of the prototype strain. This G 
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to R substitution at amino acid 766 may result in increased cleavability and the cleavage process may 


play a part in fusion activity and viral infectivity [Lau et al., 2011; Vijgen et al., 2005]. 


There were minimal amino acid changes at Spike gene nucleotide positions 448 to 459 (whole genome 
positions 24091-24102) in the glomerular part of the Spike gene, encoding the TQDG (a.a. positions 
150-153), and low mutation frequencies in that region for our four clinical isolates, although all four 
had a Y154V substitution. This is in the lectin domain of the $1 subunit involved in attachment of the 
virus to the cellular receptor which is a derivative of neuraminic acid [Kin et al., 2015]. The TQDG 
sequence may have evolutionary and functional importance, is not present in the closely related bovine 
coronaviruses and not inserted in some strains of HCoV-OC43 [Kin et al., 2015]. Four critical sugar- 
binding residues Y168, E188, W190 and H191 were present in our HCoV-OC43 isolates corresponding 


to Y162, E182, W184 and H1835 in the bovine coronavirus sequence [Peng et al., 2012]. 


The receptor-binding of SARS-CoV in the S1 carboxy-terminal domain around the receptor binding 
motif involving amino acids 479 and 487 are areas where nonsynonymous substitutions and genetic 
diversity occurred for our HCoV-OC43 isolate sequences. These areas are important for binding 
affinity to human angiotensin-converting enzyme 2 for SARS-CoV, and host immune responses [Li, 
2015]. Our isolates had nonsynonymous substitutions at amino acid residues 259 and 260-264 in the N- 
terminal domain of S1, and, in particular, had nonsynonymous amino acid substitutions within the 
region between amino acids 471 and 550 at a putative receptor binding domain of HCoV-OC43 (a.a. 
positions 339-549) [Lau et al., 2011]. In other HCoV strains, for instance HCoV-229E, antibody 
neutralization of the virus is dependent on the antigenic phenotype of the S1 subunit region [Shirato et 
al., 2012]. Hence, this may also be a factor for HCoV-OC43 needing further study. Our patients all had 
acute upper respiratory and systemic signs and symptoms. The small number evaluated here precludes 
conclusions about pathogenicity and viral mutation rate patterns, but this approach opens up pathways 


to future study. 
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In conclusion, quasi-species were present in our HCoV-OC43 strains involving areas of the S1 subunit 
of the Spike gene that may affect evolution of the viral binding process and antigenicity, as well as host 
specificity. Further studies are needed to more fully characterize the extent of quasi-species in more 
clinical isolates and their relation to disease characteristics and host immune responses. 
Nonsynonymous mutations were more frequent than synonymous ones, contributing to the hypothesis 
that the mutations are important to viral persistence in the human population over time and to disease 


pathogenesis. 
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Figure Legends 


Figure 1. Sliding window analysis of viral mutation loads, calculated as normalized Shannon 


entropy, over the entire HCoV-OC43 Spike gene (4086 bp) for the four clinical samples. 


Figure 2. Neighbor-joining trees constructed with HCoV consensus sequences. 


Panel A. All four HCoV-OC43 consensus Spike genes with seven references showing 


clustering of all four HCoV-OC43 samples with HCoV-OC43 strain (AY391777). 


Panel B. Detailed analysis of the four HCoV-OC43 samples showing two clustering in genotype 


D and two clustering in genotype B. 


Figure 3. Tight alignment of consensus amino acid sequences for the Spike genes of the four 
HCoV-OC43 clinical samples and the prototype strain, AY391777. Dots indicated the identity and 


the asterisk donated stop codons. 
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Table 1. Clinical Characteristics of Acute Respiratory Illnesses in Patients Associated with the Four Human Sequenced Coronavirus OC43 (HCoV- 
OC43) Clinical Samples 


Anti-HCoV-OC43 Serum / 
Nasal Wash Reciprocal 
Antibody Titers* 


Acute 
Acute Illness 
Study Age Date of Illness Convalescent Symptoms/Signs of Acute Severity Acute VAS of 
Subject (Years) Gender IlIness Visit Visit Respiratory Illness Score” Severity‘ Medical History 
2 39 Male 1-6-2010 664/<5 1,408/36 Spann hininisaysp nea) 45 5 Nowe 
headache, fatigue, sore 
throat 
Cough, sputum, rhinitis, COPD“, emphysema, 
3 82 Male 12-17-2010 1,983/59 6,758/1,686 dyspnea, chills, headache, 57 7 ischemic heart 
myalgia, fatigue, sore disease, diabetes 
throat mellitus 
Cough, sputum, rhinitis, Congestive heart 
d hills, headach fail ischemic heart 
4 61 Male 1-10-2011 dee [isi 80 6 Sain aass 
myalgia, body aches and disease, asbestosis, 
pains, fatigue, sore throat, diabetes mellitus 
pharyngitis 
6 62 Male 3-13-2012  1,510/139 —«:1,042/460 BRINIES yep ds Clie: 27 6 PoneMmle neat 
body aches and pains, disease, diabetes 
fatigue, sore throat mellitus, sinusitis 


“Serum antibodies were IgG and nasal wash antibodies were IgA binding to UV light and psoralen-inactivated tissue culture-adapted HCoV- 
OC43 (ATCCHVR-1558) measured by enzyme-linked immunosorbent assay. 

Severity of influenza-like symptoms and signs score. 

“VAS is visual analogue scale score. 

4COPD is chronic obstructive pulmonary disease. 
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Table 2. Numbers of Mutations in the Four Human Coronavirus OC43 (HCoV-OC43) Samples at 0.5% Mutation Frequency and 


Quality Score >30 
No. of Mutations (% of total) 
Clinical HCoV-OC43 Samples Synonymous Non-Synonymous Total 


2 9 (39%) 14 (61%) 23 
3 2 (10%) 18 (90%) 20 
4 5 (12.5%) 35 (87.5%) 40 
6 12 (32%) 26 (68%) 38 
Total 28 (23%) 93 (77%) 121 
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Table 3. Complete Listing of 121 Viral Mutations and Population Frequencies in Four Human Coronavirus OC43 Clinical Samples 


Sample 2 Sample 3 
(genotype B) (genotype D) 
Nucleotide Position in: Frequency Nucleotide Position in: 
Spike Complete Spike Complete 
Gene Genome Gene Genome 
29 23672 1.26% 41 23684 
65 23708 0.99% 65 23708 
79 23722 36.92% 81 23724 
81 23724 7.96% 174 23817 
104 23747 0.56% 259 23902 
734 24377 0.53% 372 24015 
735 24378 0.65% 997 24640 
736 24379 0.67% 1012 24655 
801 24444 0.67% 1675 25318 
894 24537 0.50% 1682 25325 
979 24622 0.97% 1691 25334 
1036 24679 0.51% 1709 25352 
1431 25074 0.52% 1712 25355 
1773 25416 0.57% 1716 25359 
2172 25815 1.05% 1730 25373 
2244 25887 2.44% 1732 25375 
2383 26026 0.51% 1736 25379 
2498 26141 0.57% 1754 25397 
2858 26501 2.10% 2000 25643 
3133 26776 0.59% 2498 26141 
3480 27123 0.50% 
3958 27601 0.96% 
3970 27613 0.68% 
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Frequency 


0.50% 
0.83% 
0.92% 
1.25% 
0.55% 
0.55% 
0.71% 
0.62% 
0.71% 
0.53% 
0.79% 
0.68% 
0.74% 
0.61% 
0.56% 
0.52% 
0.60% 
0.72% 
1.22% 
0.52% 


Sample 4 
(genotype D) 

Nucleotide Position in: 

Spike Complete 
Gene Genome 
29 23672 
41 23684 
65 23708 
81 23724 
94 23737 
98 23741 
104 23747 
119 23762 
169 23812 
242 23885 
284 23927 
310 23953 
320 23963 
325 23968 
740 24383 
784 24427 
1229 24872 
1457 25100 
1566 25209 
1581 25224 
1622 25265 
1641 25284 
1670 25313 
1675 25318 
1691 25334 
1913 25556 
2494 26137 
2738 26381 
2754 26397 
2775 26418 
2948 26591 
3079 26722 
3095 26738 
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Frequency 


1.21% 
0.84% 
0.95% 
0.82% 
0.89% 
0.56% 
0.54% 
0.51% 
0.59% 
0.55% 
0.75% 
0.53% 
0.81% 
0.59% 
0.55% 
0.58% 
2.03% 
0.51% 
0.68% 
0.54% 
0.51% 
0.60% 
0.57% 
0.62% 
0.60% 
1.07% 
0.61% 
0.50% 
1.05% 
0.61% 
0.95% 
0.77% 
0.65% 


Sample 6 
(genotype B) 

Nucleotide Position in: 

Spike Complete 
Gene Genome 
6 23649 
29 23672 
65 23708 
79 23722 
81 23724 
104 23747 
211 23854 
284 23927 
320 23963 
325 23968 
541 24184 
672 24315 
766 24409 
820 24463 
929 24572 
1217 24860 
1437 25080 
1503 25146 
1523 25166 
1544 25187 
1553 25196 
1596 25239 
1600 25243 
1601 25244 
1777 25420 
1841 25484 
1859 25502 
2078 25721 
2283 25926 
2383 26026 
2826 26469 
2862 26505 
3120 26763 


Frequency 


0.54% 
0.98% 
0.68% 
45.12% 
7.07% 
0.54% 
0.92% 
0.95% 
0.51% 
0.55% 
0.60% 
0.72% 
0.93% 
0.50% 
0.66% 
0.62% 
0.50% 
0.81% 
0.50% 
0.63% 
0.63% 
0.68% 
0.57% 
0.65% 
0.89% 
0.97% 
2.27% 
0.62% 
0.70% 
0.56% 
0.55% 
0.52% 
0.51% 


Sample 2 Sample 3 Sample 4 Sample 6 


(genotype B) (genotype D) (genotype D) (genotype B) 
Nucleotide Position in: Frequency Nucleotide Position in: Frequency Nucleotide Position in: Frequency Nucleotide Position in: Frequency 

Spike Complete Spike Complete Spike Complete Spike Complete 

Gene Genome Gene Genome Gene Genome Gene Genome 
3137 26780 0.51% 3151 26794 0.55% 
3309 26952 0.89% 3277 26920 0.80% 
3462 27105 0.66% 3284 26927 0.83% 
3676 27319 0.59% 3465 27108 0.54% 
3763 27406 1.10% 3620 27263 0.84% 
3929 27572 0.60% 
4059 27702 0.57% 
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Figure 2a 
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